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Abstract 



r> I We present numerical and analytical studies of the equation of state of translationally 

c^ I invariant black hole solutions to five- dimensional gravity coupled to a single scalar. As an 

application, we construct a family of black holes that closely mimics the equation of state of 

quantum chromodynamics at zero chemical potential. 
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1 Introduction 

In the supergravity approximation, the near-extremal D3-brane has equation of state s oc T^, 
with a constant of proportionality that is 3/4 of the free-field value for the dual A/" = 4 super- 
Yang-Mills theory [1]. The speed of sound is c^ = 1/^/S, as required by conformal invariance. 
On the other hand, the speed of sound of a thermal state in quantum chromodynamics 
(QCD) has an interesting and phenomenologically important dependence on temperature, 
with a minimum near the cross-over temperature Tc. Lattice studies of the equation of state 
are too numerous to cite comprehensively, but they include [2] (for pure glue), [3] (a review 
article), and [H [5] (recent studies with 2 + 1 fiavors). 

We would like to find a five-dimensional gravitational theory that has black hole solutions 
whose speed of sound as a function of temperature mimics that of QCD. We will not try 
to include chemical potentials or to account for chiral symmetry breaking. We will not 
try to include asymptotic freedom, but instead will limit our computation to T ^ 4Tc and 
assume conformal behavior in the extreme UV. We will not even try to give an account of 
confinement, except insofar as the steep rise in the number of degrees of freedom near the 
cross-over temperature Tc is recovered in our setup, corresponding to a minimum of Cg near Tc. 
We will not try to embed our construction in string theory, but instead adjust parameters in 
a five-dimensional gravitational action to recover approximately the dependence CgiT) found 
from the lattice. That action is 

R-l{d<pr-v{<j>) 
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We will not include higher derivative corrections, which would arise from a' and loop cor- 
rections if the theory ([T]) were embedded explicitly in string theory. 
The ansatz we will study is 

ds' = e'^{-hdt' + dx') + e^^^ . (2) 

ft 

where A, B, and h are functions of r, and (j) is also some function of r. This ansatz is dictated 
by the symmetries: we want translation invariance in the R^'^ directions parametrized by 
(t, x), and we want 5*0(3) symmetry in the x directions but not 5*0(3, 1) boost invariance — 
because boost invariance is broken by finite temperature. Assuming conformal behavior in 
the extreme UV means that we assume the geometry ([2]) is asymptotically anti-de Sitter. A 
regular horizon arises when h has a simple zero. Let's say the first such zero (that is, the 
one closest to the conformal boundary) is at r = rj^. It is assumed that A and B are finite 
and regular at r = r//. Standard manipulations lead to the following formulas for entropy 
density and temperature: 
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and once these quantities are known, the speed of sound can be read off from 
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The formula for the entropy density in ([3]) comes from the Bekenstein-Hawking result 5 = 
A/4:Gn, where A is the area of the horizon (really a volume in our case) and G^ = ^s/Stt. 
The formula for the temperature comes from Hawking's result T = /T;/27r where n is the 
surface gravity at the horizon. 

By adjusting V^cj)) one might expect to be able to recover any pre-specified CsiT), at least 
within certain limits — perhaps including that s{T)/T^ should be monotonic or some similar 
criterion (in this connection see [6])o The main aim of this paper is to characterize how 
V{(f)) translates into Cs{T) and vice versa. In section [2] we begin with the simplest possible 
case: Cs{T) constant. It translates into V{(f)) = Vqc"''^ for some Vq < and 7 related to c^. In 
section [3] we tackle the general case, exploiting a weak form of integrability of the equations 
resulting from plugging ([2]) into ([1]). In section [5] we exhibit several examples. These include 



^An earlier study [^ of thermodynamic properties of putative holographic duals to QCD starts with a 
lagrangian including an unspecified matter term. 



a particular V{(f>) whose corresponding Cs(T) curve closely mimics that of QCD. We close 
with a discussion in section El 

The results in this paper are based in large part on [8] , and aspects of them will also be 
summarized in [9]. 

2 Chamblin-Reall solutions and an adiabatic general- 
ization of them 

In a Z)-dimensional conformal field theory (meaning a CFT in D — 1 spatial dimensions plus 
one time dimension), the entropy density must obey 



s (xT 
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(5) 



simply because this expression is dimensionally correct and there is no scale other than the 
temperature that would permit a more complicated dependence. So the speed of sound is 
Cs = 1/V-D — 1. If -D > 4, then we could obtain a non-conformal theory in four dimensions by 
compactifying our CFT/) on a. D — 4-dimensional torus. (A similar idea has been considered 
in [ini [II]-) Doing so should not change the speed of sound: a planar sound wave in the 
resulting 4-dimensional theory would correspond to a planar sound wave in the original 
theory whose propagation is in the direction of the uncompactified directions. 
The A(iS'i5_|_i-Schwarzschild solution is an extremum of the action 
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and it takes the form 



where 



ds' = ^i -hde 



1- 



dx 



v-D 
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We use hats to distinguish Z) + l- dimensional quantities from 4-dimensional ones. It is easy to 
see that T oc 1/ Zh and s oc 1/ z^~^ ^ so that s oc T^~^ as the conformal field theory requires. 
Suppose we now perform the dimensional reduction described in the previous paragraph on 



the solution (jTj). In slightly more generality than we need, the Kaluza-Klein ansatz is 





df =exp{ \l^^ — ^0|>rfs2 + exp<(-J— — — — — — 0[>t/4_4, (9) 

where ds^ is a five-dimensional metric and ds']^_^ is the fiat metric on a torus T^^"^, whose 
shape we will assume to be square with side length £, so that VolT^~^ = C,^~'^. All com- 
ponents of the metric, and also 0, are assumed to depend only on the five-dimensional 
coordinates. It is assumed that £ is a constant; variation of the size of the torus is taken care 
of by the exponential prefactor multiplying ds'jj_^ in ([9]). The particular coefficients in the 
exponentials were chosen presciently to obtain a simple five-dimensional action. Comparing 
the general form (Q with the specific solution ([6]), one finds 

d.-(i^r'"'7-M*^ + <^.- + ^) e^^Q^'S^^^^^^, (10) 



where h = 1 — z^ / z^ as in ([8]). The line element flTOj) was obtained by the authors of |12j . 
but not via Kaluza-Klein reduction; instead, they considered black hole solutions to the 
equations of motion from an action like ([I]) with potentials of the form 

n0) = ^oe^^ (11) 

with Vq < 0. To see that the solutions have to come out the same in either approach, let's 
carry through the Kaluza-Klein reduction at the level of the action by plugging (jQj) into ([6]). 
After performing the trivial integral over T^~^, one obtains 



S = — 2 — / '^^^ y/~9 



i?-i(90f-V(0) 
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where V{(j)) has the form fill I) with the identifications 



n = -^^^^ 7 = a/-^^. (13) 

° L2 ^ \3D-1 

Evidently, the length scale i enters the action only as a prefactor, which can be absorbed 
into a definition of the five- dimensional gravitational constant: Kg = fi:|)^^/£^~^. 

By comparing the expression for 7 in flT^ with the result Cg = 1/V-D — 1 for the speed 



of sound, we find 

^^ = 1-2! 

^^3 2 



c: = ^--r. (14) 



This result can be derived more directly by showing that s oc T®/^^ '^"' ^ for Chamblin-Reall 
solutions: explicitly, 



. = — r-V'' = — ex 1-^ 

D 1 8-37^ r/7 1,^ 

^ = = ^ exp {[-- — ](pH 
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where (pn is the value of at the horizon. The dimensional reduction we have described 
is well-defined only for integer D > 4, but for the purposes of the computations presented 
here, it can be any real number greater than 4. 
Suppose we rewrite the result flTSl) as 

logs = h (constant in (pn) 

logT = ( 2 ~ 3~ J '^^ "^ (constant m (pn) ■ 

Given (1161) and the formula 7 = V'{(f))/V{(f)), a natural next step would be to guess the 
following dependence of s and T on (pu when 7 is a slowly varying function of rather than 
a constant: 

f^" ^(0) 
logs = - d(f) — — — + (slowly varying in (pn) 

h, ^('^) ^7N 

r't'H flV'icp) lVi(p)\ 

^°^^ ^ 1 '^'^ U vw - 3 n0) J + ('^°"^^ ""'^^"^ '" ^^^ ■ 

The lower limit cpQ in the integrals is an arbitrary cutoff. If we assume that V{(f)) has a 
maximum at = and an expansion of the form (137|) . then V{(j))/V'{(j)) ~ — 12/(?7i^L^0) 
near = 0. So the integrals in (TT7|) diverge if they are continued all the way to = 0, and 
the cutoff 00 must be chosen to have the same sign as cpu to avoid this divergence. 
A consequence of the estimates (IT7I) is a simple formula for the speed of sound: 



^2 _ d\ogT/d(PH _ 1 _ 1 V\<pH. 
^' d\ogs/d(f)H ^ S 2 V{(f)H) 
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Another consequence is 



s 3 f^" V'((j)) 

log — = -- / c/0 -— — + (slowly varying in (pn) 

3 , V{(I)h) , , , • • ^ N 

= ~ o ^°g T77TT + (slowly varying m 0^) • 
2 \/(0o) 



(19) 



A simpler way of expressing (1191) is 

" ex |y(0H)|-^/^ (20) 
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up to corrections from slowly varying terms. This is interesting because s/T is one way of 
defining the effective number of degrees of freedom available to a system, and we see from 
fl20|) that it is closely related to the potential evaluated at the horizon. 

The results ( TT8l) and ( TT9l) are a first attempt at solving the problem of translating an 
arbitrary V{(j)) to an equation of state, or an arbitrary equation of state into V{(j)). Here's 
how the latter process would work. Suppose one specifies the equation of state as s = s(T). 
Ignoring corrections to (TT9!) . one has 

/.-?log^ = log^, (21) 

where Vq is some constant. Let's regard / as the independent variable. Because V = V^e^ , 
all we need is to find = 0(/), and we will have a parametric representation of V{(j)). One 
may rewrite (IT5|) as 

""^ " 3 " 2{d<P/dff ' ^^^^ 

where corrections have again been ignored. Knowing s{T) with good precision means one 
can express c^ as a function of /. Then fl22l) can readily be integrated to give 

<i>U) = I ; • (23) 



v/2 (I - c.Uf) 



The integral is left in indefinite form because adding a constant to is obviously allowed. 

We stress that the result of plugging fl23|) into the form V = Vqc^ will result in a V{(f)) 
that only approximately reproduces the desired s{T). If the speed of sound varies rapidly 
with T, the approximation may be poor. In section H] we will show how to improve this 
approximation without resorting to differential equations that cannot be explicitly solved in 



terms of indefinite integrals. 

3 A nonlinear master equation 

There is a residual gauge freedom in the ansatz ([2]), namely reparametrization of the ra- 
dial direction. A convenient gauge choice, which should be at least piecewise valid in any 
geometry where the scalar is non-vanishing, is to set r = 0. Then the line element becomes 

ds' = e'^{-hdt' + dx') + e^^^ , (24) 

h 

and the equations of motion following from the action ([T]) take the form 

A" - A'B' + - = (25a) 

6 

h" + {4A' - B')h' = (25b) 

6A'h' + h{24A'^ - 1) + 2e^^V = (25c) 

4A' -B' + -^- V^' = , (25d) 

h h 

where primes denote d/dcj). The first two of these equations come from the tt and x^x^ 
Einstein equations; the third comes from the 00 Einstein equation; and the last comes from 
the scalar equation of motion. There is typically some redundancy in equations obtained 
from classical gravity, with or without matter. In the case of (l25l) . the redundancy is that 
the $ derivative of the third equation follows algebraically from the four equations listed. 

The ansatz flM|) has one peculiar feature: e^^ must have dimensions of length squared. 
This is because is dimensionless. 

The equations of motion ( 125|) enjoy a weak form of integrability, in the following sense: 
If a smooth "generating function" G(0) is specified, then it is possible to find a black hole 
solution where A'(0) = (7(0) in terms of indefinite integrals of simple functions of G{(j)) and 
G'{(f)). But V^(0) itself is expressed in terms of such integrals, and one cannot easily find all 
the possible G(0) that lead to a specified V{(j)). In other words, there can be simple analytic 
solutions to (123]) for special V^(0) at a special value of the temperature, but as far as we 
know, there is no nontrivial V{(j)) (i.e., none besides the exponential form) for which analytic 
solutions exist over a continuous range of temperatures. 

To understand this claim of integrability, let us consider v4'(0) = (^(0) to be fixed as a 



function of and work out A{(f)), B{(f)), h{(j)), and V{(j)). The first of these is triviah 



A{<j)) = Ao+ / rf0G(0). (26) 

J (f)0 



Computing -6(0) is immediate once one solves (I25al) for B': 



ho G{(j)} 

Next one observes that (I25bl) is straightforwardly solved once one knows A{(j)) and B{(j)): 

h{<p) = ho + hi f 4e-^^(^)+^(^) . (28) 

J<f>o 

Now (I25cp can be solved for V{(j)) in terms of known quantities: 

-2B(<t>) / h'(6)\ 

Vi^) = h{^f-^- [l - 24G(0)2 - ^^('^^T^j • (^^^ 



The constraint equation (125dp doesn't yield any new information. 
If one chooses 

G(0) = -^ , (30) 

37 

then by working through fl2Bl) - fl2^ one recovers the Chamblin-Reall solution in the form 

A(0) = Ao - 



S(0) = i?o-|(0-0o) 



37 
-0o) 

37^ 



(31) 



where 



/i(0) = ho + hi exp 

I 07 

y(0) = Voe^* 



h = , ,J /^i Vo = -^-^e-^^o--^o/.o . (32) 

8 — 37^ 67 



By choosing 

one obtains V{(j)) in a form that doesn't depend on any integration constants at all. This 



0o = -(log/io-25o) , (33) 
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situation is very special and corresponds to the fact that for V{(f)) oc —e'^'^ one can find a 
whole family of analytic solutions parametrized by (pH- 

By differentiating combinations of fl26|) - fl29l) . one can derive the following "nonlinear 
master equation:" 

G' d /G' 1 ^^ G' \ .^^. 

Xllog T7 + 7?7^-4G-— — — - . (34) 



G + V/3V' d(f) ^\G 6G G + V/3V' 

Describing ( 134|) as the master equation is appropriate because if one starts knowing V{(f)) 
and manages to solve flMl) . then to obtain a black hole solution one only needs to perform the 
"trivial" integrations (!26l) - (!28l) . A numerically efficient strategy for obtaining an equation of 
state given V{(j)) centers around solving (134|) numerically. In more detail, the procedure is: 

1. Choose the value (pn of the scalar at the horizon. 

2. Find a series solution of the nonlinear master equation around (p = (pu- 

3. Seed a numerical integrator like Mathematica's NDSolve close to = cpn using the 
series solution. 

4. Integrate the nonlinear master equation up to a value of cp close to a maximum, corre- 
sponding to the asymptotically anti-de Sitter part of the geometry. 

5. Extract s and T from integrals of simple functions of G{(p). 

Most of these steps require further explanation, which will occupy the remainder of this 
section. 

At the horizon, h has a simple zero, and the other quantities appearing in fl25cl) and 
(125dl) are finite. Evaluating these two equations at the horizon gives 

V{<Ph) = -3e-2^(^-)G(0H)/i'(0H) V'{<Ph) = e-2^(^-)/i'(0^) , (35) 

which implies that G + V/3V' vanishes at the horizon. Starting from this condition, one may 
develop a power series solution around the horizon: 

This series solution can be developed to any desired order without encountering arbitrary 
integration constants. 



To understand the asymptotic behavior far from the horizon, let's speciahze to the case 
where V{(j)) has a maximum at = 0: 

12 1 
V{<P) = -- + -m'<P' + 0{<P'), (37) 

where m'^ < in order for = to be a maximum. The gauge theory operator O^ dual to 
(f) has dimension A, where 

A(A - 4) = m^L^ . (38) 

We will be primarily interested in the case where A is close to 4. It helps our intuition at 
this point to pass to a more standard gauge: instead of setting r = 0, we can set i? = to 
obtain 

ds^ = ^^{-hdt' + dx'') + '^. (39) 

We note that A and h appearing in fl39|) are precisely the same as when we use the r = (p 
gauge, only expressed as functions of r rather than 0. Large r corresponds to the region far 
from the horizon, and the leading asymptotic behavior of solutions there is 

A^y h^l ^ (AL)^-^e(^-^)^ . (40) 

Each approximate equality in (140|) means that the ratio of the two expressions on each side 
approaches 1 as r becomes large. The behavior of indicates a relevant deformation of the 
conformal field theory: 

£ = /:cFT + A'-^0^. (41) 

As we vary temperature to compute the equation of state, we should of course keep A fixed. 
A simple way to do this is to set AL = 1. Then the last equation in (140!) becomes 

^(0) = ^^ + o(l) for small 0, (42) 

where o(l) means a quantity that is parametrically smaller than 1 in the limit under 
consideration — so in the limit — ;> in the case of (142|) . 

In order to compute the entropy density using Q, we need to know A(0//). This can be 
extracted by comparing (H2I) to (I2BI) with 0o set equal to (f>H and Aq set equal to A^ = A{(J)h)'- 

A(0) = Ah+ f d^G{4>) = ^ + 0(1) . (43) 

J<i>H ZA-4 
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Solving for Ah and then taking — i> 0, one finds 



A 



H 



log</'H 

A-4 



G(0)- 



(A-4)0, 



(44) 



Tlie integral converges because the explicit 1/0 term cancels against the leading behavior of 
G(0) for small 0. Plugging (1H|) into the expression for entropy density from ([3]), we find at 
last 



27r 



3/(A-4) 



^2'1'H 



exp < 3 



G(0) 



(A -4)0, 



(45) 



A similar formula for the temperature may be derived starting with the observation that 

^ = -^^ (46) 



where B is the function controlling the g^^i, metric component in r = gauge. One obtains 
fH6|) by comparing the flM|) to fl39|) . The sign is based on assuming that increases from 
to positive values as r decreases from +oo to finite values. The first equation in fHUl) implies 
dA/dr — ^ 1/L as r — i> oo. Combining this with ( l46l) gives 



G 



dA dr dA 
dd> dd) dr 



R 



L' 



(47) 



where the approximate equality means that the ratio of the last expression to the previous 
ones approaches 1 as goes to 0. In summary, 



1 ^ -LG(0)e 



-BW 



(48) 



using the same sense of approximate equalities. (Recall that e~^ has dimensions of inverse 
length, while G(0) is dimensionless, so (H8|) is dimensionally correct.) We may cast the ex- 
pression for temperature from ([3]) in terms of G(0) by multiplying by 1 in the form indicated 



m 



^ ^ e^«-^(M\h'{^H)\ ^ e-^^-^(^^)/^-(0^) ^ .(,, 

Le-^B(<t>H)G{(j)H)h'{(t)H 



An 



exp <^ Ah + B{<J)h) - B{(j)) - log 



g(0f/) 

G(0) 



(49) 



LV{(t) 



H) 



1271 



exp <Ah + 



rf0 1 
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In the second step we assumed that h'{(j)H) < 0, which is the expected sign when vanishes 
on the boundary and is positive at the horizon. In the last step we used the first equation 
from fl35l) to simplify the prefactor and also 

B(M-BW = lo,^+r-^. (50) 

G(0) J^ QG{(j)) 

which is a consequence of ( 1271) . The integral in the last expression of ( l5T1) converges even 
when — i> 0. We use (1H|) to eliminate A^ from (15T|) and obtain at last 



i = —^ -^exp 



Gi4>) - 77 JTT + 



1 



(A -4)0 6G(0)_ 
The measure of the number of degrees of freedom that is easiest for us to access is 



(51) 



where to obtain the right hand side we simply combined (l45l) and (ISTj) . When (pH is small, 
entropy and temperature become large because of the factors (pjj ~ and (f)j^ ~ in (l45l) 
and (151 p . In this limit, the integrals in the exponent become negligible, and (1521) becomes 



•S 


4 


p 



(53) 
So we recover conformal behavior in the ultraviolet, as expected. 

4 An approximate determination of the equation of 
state 

The adiabatic formulas (IT7|) work well when (pn is in a region where V{(j)) is nearly exponen- 
tial, but they do not work well for small (pn, where ^(0) is close to attaining a maximum. 
This is shown in figure [T] for V{(f)) = — (12/L^) cosh(0/2). On the other hand, it's easy 
to extract asymptotic formulas valid in the (pn —^ limit from (HSj) and (IFIl) : using the 
expansion (!37|) . one finds 
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d(log T)/d<^^ 






12 
V(^)= — rcosh- 



d(log s)/d(^^ 
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Figure 1: A comparison of the exact ^^^^ and "^^^"^ '' -' for V^(0) — -p^ ^^^,, ^ 
adiabatic approximation, flTTl) . and the improved approximation scheme, fl57|) with the choice 



We wish to find formulas that interpolate smoothly between (TT7j) and (!Ml) and involve at most 
indefinite integrals, not solutions to a difficult, nonlinear, second-order differential equation 
such as ( 1M|) . We start by noting that the formulas fl5^ together with fl37j) imply 



(55) 



(ilogT A/ir(0H) 1^0^) 



o?logs _ A /^ V{(^h) 
~ 4 



^/f 



y' 



'H, 



for small </)i:f. 



where approximate equality means that the ratio of the two sides tends to 1 as 0^ -^ 0. On 
the other hand, provided V^(0) tends to an exponential form V^e^'^ for large 0, the adiabatic 
approximation becomes good if ^n is large. So for such potentials, (ITTIl can be rephrased as 



d<PH ^ 2 V{<Ph) 
dlogs _ V^(0g) 



V 



for large 0^- 



(56) 



>H 



Comparing (l55l) and (156|) . we are led to introduce "fudge factors" Ps{(Ph) and pt{4>h) such 

that 

dlogT ^ ,^ ^ A^'(</'/^) 1^(0/^) 



acpH 
dlogs 

d(t)H 



Pt[(Ph) - 



Ps[(Ph) 



2Vi(l)H) 3\/'(</.^) 
V{<I>h) 



(57) 
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We can rephrase (155!) and (156!) as the statements that both Pt{<Ph) and Ps{4>h) interpolate 
between A/4 at small (pn and 1 at large (pn- Our improved estimate of the equation of state 
consists simply of guessing an interpolating function with these properties. The guess is 

PTi<l>H) ^ Ps{4>h) ^ pi<l>H) = 1 + ^^ (t - ^) • (^^) 

In (1581) approximate equality means that pT, Ps, and p are supposed to be nearly equal for all 
(pH- But away from the small (pn and large (pn limits, fISS]) is not a controlled approximation, 
in the sense that there isn't a parameter that we can tune to make it better. It is nevertheless 
useful for quickly determining the qualitative features of an equation of state given V{(p), 
as illustrated in figure [H There might be a better choice of Pt{(Ph) and Ps{(Ph), or perhaps 
even a systematic expansion for them in terms of powers of the potential and its derivatives. 
Starting from (IFr|) . we have immediately 

2 d\ogT/d(PH Pt{<Ph) fl IV\(Ph) 



^ d\ogs/d(PH Ps{cPh) V3 2V(0^)2; • ^ ^ 

Thus, assuming Pt{4>h) ~ Ps{(Ph) is the same as assuming that the speed of sound, as a 
function of (pn, is well-approximated by the adiabatic formula, flTSl) . 



5 Examples 

The simplest analytical form that interpolates between a maximum a.t (p = and exponential 
behavior for large (p is 

V{(P) = V,oM) = -^cosh7(/.. (60) 

The adiabatic treatment discussed in section [2] leads us to expect that the speed of sound will 



he Cs = -i/| — ^ for large (pn- So in order to have stable black holes in this regime, we must 
have 7 < y2/3. This bound can be regarded as an application of the correlated stability 
conjecture (CSC) [131 [T3]. But there is a tighter bound on 7 coming from the behavior near 
= 0: 

Kosh(0) = -^- ^0^ + 0(0^), (61) 

so m^ = — 127^/L^. In order to satisfy the Breitenlohner-Freedman bound, m?L'^ > —4 
p!5| [T6| [T7] . we must have 7 < 1/V3. This is very restrictive, because it means that the 
minimum speed of sound we can arrange at large (pn using the pure cosh potential (J60l) is 
Cs = yS/lS ~ 0.53, only slightly smaller than the conformal value Cs = l/v3 ^ 0.58. The 
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Figure 2: The speed of sound for V{(f)) = — jf cosh -4. 

behavior of c^ as a function of T is shown in figure [2] for 7 = 1/Vq. 

If one uses the potential (pU]) . then Cs in the infrared is tied to the dimension A of the 
operator that breaks conformal invariance in the ultraviolet. Let us consider a minimal 
generalization that loosens this artificial constraint: 



v{4>) = v{j,b-. 



12 
"Z2 



cosh •ycf) + 6</)^ 



(62) 



A parameter equivalent to 7, as before, is the speed of sound in the infrared, c^ = | — -2"- 
With 7 fixed, a parameter equivalent to b is the dimension A of the operator dual to </> near 
the UV fixed point: 



Note that taking A close to 4 amounts to making V{(f)) nearly quartic near its maximum. 
As we will report in more detail in [9], the choice 



7 = 0.606 



2.06 



(64) 



corresponding to c^ = 0.15 in the infrared and A = 3.93, leads to an equation of state 
that bears a striking resemblance to the one expected for QCD: see figure |3l It may seem 
surprising that the most distinctive feature of the equation of state of QCD, namely a 
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lattice data vs. black holes 
0.35 f 2+lAPPROX pure glue 



cl 




Figure 3: The equation of state of a black hole (red) compared to the lattice equation of 
state for pure glue (blue) and 2 + 1 QCD. The pure glue curve is based on [2] and private 
communications from F. Karsch. The 2 + 1 QCD points are based on [5]. 

smooth but rapid cross-over, emerges from a potential that is nearly featureless. To gain 
some intuition about why this happened, consider again the adiabatic approximation (ITSl) 
to the speed of sound. When (pn is close to where the nearly quartic behavior of V{(f)) rolls 
over into nearly exponential behavior, this approximate formula predicts that c^ dips to a 
fairly low value, only to rise back up again for larger (pn towards its infrared limit, 0.15. See 
figure m 

Other behaviors emerge from the potential (!62|) for other choices of b and 7. For example, 
if 7 > y2/3, the adiabatic approximation suggests that there is a minimum temperature 
Tmin for black hole solutions. A particular case is illustrated in figure [5l Solutions with 
very low entropy have high temperature and negative specific heat, and they are always 
thermodynamically disfavored compared to a branch of high-entropy solutions. Presumably 
there is a first order transition to geometries with no horizon, similar to the Hawking-Page 
transition [18]. This transition probably happens at a temperature above Tmin. It is worth 
noting that the specific heat diverges at T^[^ because T reaches a minimum as a function of 
(pH while S is varying smoothly. 

It is also possible to have a first order transition between high entropy and low entropy 
black holes. An example where this happens is illustrated in figure [61 For a finite range 
of (pH, the speed of sound is imaginary, indicating a Gregory-Laflamme instability. This 
touches once again on the correlated stability conjecture (CSC), so let us pause to review 
it. It was proposed in [131 [E] and further argued in [19] that, in the absence of conserved 
charges related to gauge symmetries, existence of a Gregory-Lafiamme instability [20| [21] 
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Figure 4: Left: The potential (162|) with the parameter choices ( !63|) that give an equation 
of state resembhng QCD's. Right: Although V^(</)) is relatively featureless, the adiabatic 



formula c? 



1 iv'{<t>Hr 



3 2 V(^Hf 

fc 



suggests that the equation of state resulting from it will indeed 



exhibit a low minimum for the speed of sound. 
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Figure 5: The equation of state for V{(j)) = — -p- cosh \/j(p + ji'P'^- 
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is equivalent to positivity of the specific heat, C = TdS/dT. According to a more general 
version of the CSC, dynamical stability of a horizon is equivalent to positivity of an appro- 
priate Hessian matrix of susceptibilities, one of which is the specific heat [131 [H]. As pointed 
out in [22], C > is equivalent to c^ > 0, which makes the CSC seem inevitable, at least 
in the absence of conserved charges and in the presence of some kind of holographic dual. 
The argument of [22] can probably be extended to cover the general case by considering the 
dispersion relations for all the hydrodynamical modes, including those arising from the dual 
conserved currents. However, the CSC remains a conjecture, and there appears to be room 
for violations: see for example [231 IMl 122] • 

The CSC relates only to the existence of a linearized instability around a static or sta- 
tionary horizon. Considerable work has been devoted to the question of what the endpoint 
of the evolution of the Gregory-Lafiamme instability might be: see for example [2^1 [271 ESI 
|221[3niEI]- When there are thermodynamically stable horizons both with larger and smaller 
entropy, it seems to us likely that the endpoint of the evolution is a mixed phase with uni- 
formly small curvatures outside the horizon, which remains unbroken. A mixed phase is a 
configuration where high entropy and low entropy regions with the same temperature are 
separated by domain walls. Typical solutions may not be static, but may instead evolve 
slowly according toward larger domains according to an effective theory with domain walls 
whose width is eventually negligible compared to the size of the domains. Mixed phases were 
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Figure 7: The equation of state and the speed of sound for V^(0) = — |t cosh -j= + ^^0^. 
The point where c^ = is a second order phase transition. If one considers instead V{(j)) = 
— -p cosh -j= + h(t)^ for h > 1.942, the transition becomes first order, while if 6 < 1.942, it is 
a cross-over. 

previously suggested in connection with the Gregory-Laflamme instability in [32] . 

Finally, it is possible to arrange second order behavior by tuning the potential V{(f)) so 

that cj. goes to at some value of (pH but never becomes negative: see figure [71 There is 

a corresponding critical temperature, and near it the equation of state typically takes the 

form 

T — T 

s ^ sq + Si/at^/^ where t = — — — ^ . (65) 



The specific heat diverges as C ~ t ^/^, and consequently the speed of sound behaves as 

|f|-l/3 

6 Discussion 

Since the inception of the anti-de Sitter / conformal field theory correspondence [331I3HI35] . 
it has been hoped that it would help solve quantum chromodynamics (QCD). This hope was 
articulated most clearly in the early literature in [36]. Subsequently, a large and somewhat 
heterogeneous literature has grown up around the idea of "AdS/QCD." Points of entry into 
this literature include [371 [38l [391110] . 

The first thermodynamic question one might ask of a putative dual to QCD is whether 
the equation of state is right. We have shown that the equation of state can be built into the 
construction by choosing an appropriate potential V{(f)) for a scalar field that describes the 
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breaking of conformal invariance. Indeed, within certain limitations, any equation of state 
s = s(T) can be translated into a choice of V{(f)), and vice versa. The limitations include 
that we use the supergravity approximation. This immediately points to a weakness of our 
approach: the shear viscosity will always satisfy rj/s = l/47r, regardless of temperature [HI 
H2I B3l m] . Low shear viscosity is in conflict with expectations for the low-temperature phase 
of QCD, where the mean free path becomes large. Another reason to be suspicious of any 
attempt to describe the low-temperature phase using a black hole horizon is that at large N, 
entropy of a horizon scales as A^^, whereas the number of degrees of freedom in the confined 
phase of an SU{N) gauge theory scales as A^°. A black hole description may be approximately 
valid above T^, and its validity may fail only gradually as one passes through the cross-over. 
But sufficiently far below the transition, the paradigm of weakly interacting hadrons should 
take over, and that is not part of our construction. One might imagine improvements on 
our construction, where, for example, higher curvature corrections significantly increase rj/s, 
especially around or below the transition temperature. Eventually — perhaps when curvatures 
near the horizon become sufficiently large compared to the string scale — there could be a 
cross-over to a gas of strings in a curved spacetime. 

Our methods for constructing black holes are more general than the particular problem 
of mimicking the equation of state of QCD. Smooth cross-overs, second-order transitions, 
first-order transitions, and perhaps even mixed phases may all be accommodated within the 
framework we have proposed. Our nonlinear master equation approach is special to the case 
of a single scalar, and it takes advantage of a weak form of integrability of the underlying 
equations. However, it is straightforward in principle to work with multiple scalars as well 
as with gauge fields: in this connection see for example ^51 HGj. It seems likely that black 
holes in suitably designed theories exhibit a remarkable diversity of phase transitions. 
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